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Abstract 



Oh| We study the main properties of the one-loop vacuum polarization 

(— I ' function ijl^u) for massless spinor QED/^ in a slab, namely, with fields 

defined on M C M^+^, such that M = {(rro, X3)|0 < 3:3 < e}, 
and bag-like boundary conditions on the boundary. We evaluate the 
J> I induced charge density and current due to an external constant electric 

field normal to the boundary; we also study the effective action for a 



, purely transverse field, identifying its e-dependent contribution. 

T— I ■ 

ff^ ■ In the presence of nontrivial boundary conditions, Quantum Field Theory 

^ . models may give rise to many interesting effects. Noteworthy examples of 

them are the Casimir effect the bag model of QCD as well as many 

others [1]. 

The common origin of those effects is the presence of boundaries, which 
strongly affect the structure of the vacuum fluctuations; this fact has rel- 
5^ , evance not only for global observables, i.e., Casimir energies, but also for 

local ones, like vacuum energy densities and response functions. The latter, 
which are determined by the correlation between fluctuations should exhibit 
a strong dependence with the distance to the boundary (at least for the ones 
involving degrees of freedom affected by the boundary conditions). 

In this letter, we consider the vacuum polarization tensor, II^j,, for a Dirac 
field confined to a slab-shaped region A^, defined by the condition < X3 < e. 
More concretely, we consider a massless Dirac field in 3-1-1 dimensions, which 
satisfies 'bag' (i.e., vanishing normal current) boundary conditions on the two 
static planes 0:3 = and X3 = e. Euclidean spacetime coordinates are, in the 
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conventions that we shall use, denoted by (xq, Xi, 0:2, X3). Moreover, we shall 
assume the Dirac field to be confined to the region between those mirrors; 
hence, 11^^ vanishes identically outside of < X3 < e. It there are fermions 
outside of the slab, 11^^ does not vanish there. However, their properties in 
that region are equivalent to the case of just one boundary; for example, if 
X3 < 0, U^i, in that region is the same as the one for a single wall at X3 = 0. 
In the conventions we shall use, Euclidean coordinates are denoted by x^, 
/i = 0, 1, 2, 3, while the metric tensor is given by Qfj^u = Sfj,^- 

To account for the effect of the confined fermion fiuctuations on the gauge 
field propagation, one introduces: 

where the fermionic action, 5*/ accounts for the (minimal) coupling to the 
gauge field as well as for the introduction of the bag boundary conditions. 
To deal with the latter, we follow the approach of representing them by local 
interaction terms [5]. Of course, the resulting propagator should agree with 
the one one would get by using, for example, the multiple refiection expansion 
(MRE) O [71 [8]. That approach has been used in [9] for the calculation of 
U^i, for a Dirac field in a half-space. 

Following [3], we include a 'potential' V into the fermionic action, so that: 

5/(^,^;A) = J d^xilj{x)[^ + gVix3)+te^{x)]^{x) (2) 

with 

Vixs) = [5{xs) + 6{xs - e)] . (3) 

Here, 5^ is a constant which, in order to enforce bag boundary conditions, 
has to be equal to 2 0. The slash notation denotes contraction with the Eu- 
clidean 7-matrices which, in our conventions, are all Hermitian, and satisfy: 
{7m , 7m} = 2 yu, = 0, 1, 2, 3. 
From ([1]), we may write, 

Tf{A) = -Ti\n[l + iei^ + 2V)-^ 4] , (4) 

which allows us to introduce the vacuum polarization, H^u, the kernel de- 
termining the form of the first non-trivial term for the expansion of in 
powers of A: 

Tf{A) = ^ fd^x f dSA^{x)U^,{x,y)A,{y) + ... (5) 



^See [5]. Different values would produce 'imperfect' boundary conditions 
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From we may give a more explicit, yet formal, expression for 11^,^ in 
coordinate space: 

n^v(a;,y) = -e^tr [S'f(?/,x)7^5'f(x,2/)7^] , (6) 

where Sp denotes the (coordinate space) exact fermion propagator for e = 0, 
namely, the free propagator with bag boundary conditions. This object may 
be obtained as the matrix elements of the inverse of an operator: 

SAx.y) = {x\{^ + 2V)-'\y) , (7) 

where a 'bracket' notation has been used to denote the coordinate space ver- 
sions of operators (in this case, the inverse of ^ + 2V). Since the boundary 
conditions only affect the X3 coordinate, there is symmetry under transla- 
tions in the 'parallel' coordinates x„ = {xq,Xi,X2)- Thus, Sp = Sf{x„ — 
I/,,; X3, 1/3), and it is natural to work with Sp, a mixed Fourier representation 
of Sf which results by transforming the parallel coordinates only, so that 
Spixu - y„; X3, ys) SpiPu] X3, ys). 

This reduces the problem to a one-dimensional one, where Sp satisfies 
the equation: 

[73 - ^ + 2V'(x3) S^f(p„; 3:3, 1/3) = 5(x3 - 7/3) . (8) 

The solution to the inhomogeneous linear equation above may be written as 
the sum of two terms: 

SpiPu^Xs^ys) = SP{p,;x3,y3) + U{p„;x3,y3) , (9) 

where the first of them, 5*^'' , is the free fermion propagator in the absence of 
boundaries, while the other, denoted by U accounts for the boundary (bag) 
conditions. 

The explicit form of Sp^ may be found quite straightforwardly: 

S^F\p-,X3,y3) = ^[sgn(x3-y3)73-^^] e-IP"ll"l . (10) 

Since Sp^ is the Green's function of the free Dirac operator we derive, 
from dS]), an equation for U{p,„X3,y3): 

(73^x3 + « ^11) f^fen 3^3, 2/3) = -'2V{x3) [SP{p„;x3,y3) + Uip„;x3,y3)] , 

(11) 
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which imphes: 



-2 



(0), 



+ S]^'{p,;x3,e) SP{p„;e,y3) + S)p!'{p„;x3,e)U{p„,e,y3) (12) 



(0), 



One can then obtain U(p„; X3,y3) from the equation above, for example by 
first finding U (p,,, 0, y^) and U (p,,, e, y^): these can be found by evaluating ( TT^ 
at = and X3 = e, respectively, and then solving the resulting system of 
equations. 

Once those objects are found, the outcome of this procedure is an ex- 
plicit expression for U{p„; x^^ys), which can be expressed as a sum of terms, 
distinguished according to its 7-matrix content: 

U{Pu;x3,y3) = Uo{p„;x3,y3)I + Ui{p,;x3,y3){-i-—) 

\P"\ 

+ f^2(p„; 3:3, 2/3) ( - ^ j^^s) + f^3(p,„ 3:3, 2/3)73 , (13) 

where we have introduced four functions (a = 0, 1, 2, 3). The explicit form 
of these functions, for < X3 < e and < 2/3 < e, is the following: 



f^o(Pn;a;3,2/3) 



2(e2bi.k + 1) 



^v«\{^i+yi) _|_ g-biil(a;3+2/3-2€) 



(14) 



t^i(p,i;a;3,2/3) 



20 



,2|p„|e 



1) 



glPiiK^^S-J/s) _|_ ^-\p«\ixo,~yo,) 



(15) 



^2(p„;a;3,2/3) 



2(e2|Pnk + 1) 



gbii|(a;3+?/3) _ g-b„|(x3+?;3-2e) 



(16) 



f^3(p„;a:3,2/3) 



2(e2|P"k + 1) 



g|Pii|(2:3-y3) _ g-|Pii|(a:3-y3) 



(17) 



A lengthy but otherwise straightforward calculation shows that the bag 
boundary conditions are fulfilled, namely, the following equations are sat- 
isfied: 



lirn (/ - 73)S'ir(p„;x3,?/3) = 

3:3->-0+ 

lim (/ + 73)5^F(p„;a:3,2/3) = 

X3-5>e- 



;i8) 
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(as well as the ones corresponding to approaching the boundaries with the 
^3 coordinate). 

We now analyze the UV properties of the propagator Sp, as presented 
in (Q. We see that the two terms have a quite different behaviour. Indeed, 
the first one, Sp\ being the free propagator in the absence of boundaries, 
does have the well-known UV behaviour (~ in momentum space). In the 
'mixed' Fourier representation, where it depends on p„ and X3, y^, that be- 
haviour translates into a ~ \p„\^ behaviour when = y^, and an exponential 
decay when X3 7^ y^. The U term, on the other hand, decays exponentially 
everywhere, except when both X3 and y^ approach one (the same) boundary, 
namely, when either X3 -|- 1/3 — > or X3 + ?/3 — >■ 2e. Moreover, we can trace 
the origin of that behaviour to see that it comes from the Uq and U2 terms 
{Ui and t/3 always decrease exponentially, regardless of the values of X3 and 

ys)- 

As a consistency check for the expression of Sp , one can find an alterna- 
tive representation, obtained by expressing appears as 
a common factor in the expressions for Ua, as a series: 



^ 00 
n=0 



-2(n+l)|p„|e 



(19) 



As a result, and after some algebra, we obtain: 



SF{P,]X3,y3) 



+00 



SPip,; 2^3,2/3 + 2ne) 



+ 73 SP{p„; X3, -2/3 + 2ne) 



(20) 



which may naturally be thought of as an MRE representation of Sp- It is 
possible to pinpoint also here the part of the propagator that control its UV 
behaviour. Since all the terms in the sum can be written as functions of Sp\ 
it is quite straightforward to see that the terms that control the large- Ip,,! 
regime are: n = (for X3 = ys or X3 + ys = 0) and n = 1 (for X3 + ys = 2e). 

To obtain a more explicit expression for 11^^, we insert the results derived 
previously for Sp into We first note that, since 11^^ shall also be a 

function of (x,, — y„;x3,y3), we may write its Fourier transform with respect 
to the parallel variables, as follows: 



— e 



d^p„ 
(2^ 



tr 



^^(Pii; 2/3, X3)-ff,SF{Pu + K] X3, 2/3)71. 



:2ii 



5 



Thus, 



LUi 



tUL 



UUi 



(22) 



where: 
nJ,^(/c„;x3,Z/3) = 

n^,^(/c„;x3,y3) 



li'i^ik.-x^.y^) 



tr 



(27r)3 
(27r)3 
(27r)3 



SPiPu;y3,X3)'yf,S'p!'{p„ + K;x3,y3)'y^ 



tr 



tr 



SF\pu;y3,X3hij,U{p„ + K;X3,y3)'J^ 



U{pu;y3,x3ht,sP{p„ + /c„;x3, 1/3)7^ 



(23) 



(24) 



(25) 



(2vr)^ 



tr 



U{pu;y3,X3)'j^U{p„ + k„;x3,y3)'ju ■ (26) 



Inserting the exphcit expressions for 5"^'' and U presented in the Appendix, 
one sees, after evaluating the traces, that 11^^ + iP^ vanishes identically. 
Thus, 

U^,{K;x3,y3) = nJ,(fc„;x3,i/3) + n;^,(/c„;x3,z/3) , (27) 

where: H^, = fi^,^ and fijf, = n^f^. 

Before evaluating the above for some particular cases, we calculate a 
magnitude corresponding to a related effect: the boundary conditions at 
and e do break the chiral symmetry of the (massless) unconfined theory. A 
quantitative and local measure of that violation, based on the expectation 
value of a bilinear observable is the fermion condensate p{x) = {ip{x)ip{x)). 

In terms of the fermion propagator, we see that p may be expressed as 
follows: 

p(x) = -ti[SF{x,x)] . (28) 

Furthermore, taking into account translation invariance in the parallel coor- 
dinates, and the specific form of Sp'- 



p{x) = p{x3) = - 

which may be exactly calculated: 



{2ny 



tr [S'f(Pii;x3,X3 



(29) 



, , TT^ 3 + COs(^' 

P(^3) = -7T-T 



2e3 sin^i^ 



(30) 
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Figure 1: e^p as a function of x^/e. 



which is of course symmetric with respect to X3 = |, and finite everywhere, 
except at X3 = 0, e. 

As it can be inferred from Figure 1, for finite values of e, ^(xs) diverges 
on the boundaries, and has a maximum of — vr^/e'^ when X3 = e/2. It can also 
be shown that it tends to when e — )■ 00 and X3 is far from the borders. This 
concentration around the boundaries does also show up in the calculation of 
the induced vacuum current, j^, resulting from some special electromagnetic 
field configurations. Those induced currents are, in this linear response ap- 
proximation, determined by 11^^ and the gauge field corresponding to the 
given electromagnetic field. In our conventions: 

= e (^(x)7^^(x)) 

= j dSll^,{x,y) A,{y) , (31) 
and in the Fourier representation: 

UKM = -,/<i,3n„„(fc„;X3.y3)I„(*„;!/,). (32) 

We have evaluated the current for the particular case of a static electric field 
with normal incidence (xs direction). This corresponds to Fqi = where 
is a constant. In the ^3 = gauge: Aq = —Ex^, and = for /i 7^ 0. 
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Thus, 

jf,=jf,{x3) = j dy3^f,o{0;x3,y3)y3 , (33) 
and its usually more convenient real-time (Minkowski) version, j^^, becomes: 



J^(^3) 



E / c/?/3noo(0;x3,?/3)?/3 



(34) 



and, (with / ^ 0): 



JMix3) 



E / rf2/3nio(0;X3,?/3)y3 



(35) 



Since the system is parity conserving, ji = j2 = for this external field. 
Thus, we concentrate on and jfj. Recalling (127|) . we see that these cur- 
rents receive two contributions: 



JMix3) 



E / dys noo(0;x3, 1/3)1/3 



ys 



(36) 



and analogously for jf j. 

To otain j]^, we note first that: 



noo(0;a;3,?/3) 



(2vr) 



d'p e2HI-3-y3l2(i_4 

p2 



67r2|x3 -?/3|3 ' 

(37) 

which is finite for all ^ y^. A divergence appear when integrating over 
1/3, to calculate the current. However, the usual renormalization conditions 
imply that one has to use Hadamard's finite part in that integral. 
On the other hand, for n[j^Q we find the expression: 



n^o(0;a;3,2/3) =-Mi 



X3 - ys 



X3 + y3 



(38) 



where 



Miiu) 
M2(v) 



1 , C(3) 



72 



72 87r2 
+ E2(v) 



(39) 
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with 

Fiiu) 



1 



cos(7ru) 



(1 _^,)^(2) (^i 
1 



3847r2 



(1 - u)i)^ 



(2) i _ ^ 



2. 



48 



■ CSC 7rt> 



2sin(27rt;) + 7r(l 



(2) 



U 



1 + 

3 + cos(27rt> 



u 



(40) 



where C, is Riemann's zeta function and 

In Figure 2, we plot the profile of — as a function of — that follows 

from the previous results. 

For under the same external field configuration, we have found a van- 
ishing result, namely, = 0. It should not be surprising that the normal 
current vanishes at the borders, since the bag conditions should do pre- 
cisely that. The vanishing of the current inside the slab, however, should 
be regarded as a steady state configuration feature, reached after the charge 
density has adopted the profile above. 

The profile for the charge density above corresponds to the case of a 
normal electric field; to explore a situation that is, somehow, opposite to 
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that one, we presenting here an exact expression for the effective action, Tj, 
in the case where the gauge field depends only on the parallel coordinates, and 
moreover that field is a function of only the parallel coordinates. Following 
the approach of j5] for the case at hand, an entirely analogous procedure to 
the one used there, yields now: 



-r/(A||) 



det /C 



(41) 



where K. is the object 
/C = 

with 



n 



l + V {V- 73)e"^^ 
{V + 73)e"^^ 1 + V 



V 



u 



(42) 



(43) 



Here, p\\ is the 'parallel Dirac operator': p\\ = '~faDa, ol = 0,1,2, and 
D = + A||(x||). 

Using some algebra, the determinant above may also be written in the 
equivalent form: 



det/C = [det{l + V)Y det 



1 + 73 / 
1 ^ 1 



V)e 



-2tH 



(44) 



As a consequence, rj(A||) receives two contributions, one coming from the 
determinant of {l + V), which is essentially two dimensional and independent 
of e, plus an extra term F^, which does depend on e, it is finite (because of 
the exponential factor), and tends to zero when e — )■ oo: 



r.(A|i) 



-Tr log 



i-l±^n 



V)e 



-2en 



(45) 



We end this note by presenting our conclusions: We have found the form 
of of the induced vacuum current due to a normal electric field, which shows 
that the induced charge density distributes itself in order to counterbalance 
the external electric field. We obtained an expression for the effective action 
that follows from considering a parallel gauge field configuration. It contains 
a term which represents the finite-width contribution to the effective action. 
It is finite, and it contains a ^-^^ factor, which is a projector which accounts 
for the suppression of part of the fermion modes because of the boundary 
conditions. 
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